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Abstract 

We theoretically derive the Hamiltonian dynamics of dispersion-managed 
coupled breathers in optical transmission system as a result of phases-plane 
dynamics of the system. We analytically show that the contribution of a per- 
turbation term caused by a noise sourced by amplifiers in optical fiber can 
influence the amplitudes and chirps of the pulses dynamics. We also find that 
the coupled breathers dynamics depent on a certain relationship of the chirps 
and amplitudes. The results of the perturbed coupled NLS are then reduced 
to that in the unperturbed system and the single NLS case. 



* hendry_elim@hotmail.com 
t freddyzen@hotmail.com 



1 



I. Introduction 



A dispersion-managed (DM) optical fiber is designed to create a low (or 

even zero) path-averaged dispersion by periodically alternating dispersion sign 
along an optical line which dramatically reduces pulse broadening. Recently 
dispersion management has become an essential technology for development of 
ultrafast high-bit rate optical communication lines. There arc some effective 
methods to improve the performance of a soliton transmission system such as 
the use of a dispersion-compensation technique, the use of dispersion-managed 
technique and so on. 

Analytic techniques about the dispersion managed system of the single NLS 
by Kutz et. al.^and Gabitov et al.^ have been separately developed to describe 
the amplitude and phase fluctuations that occur on the scale of the dispersion 
map. The aim of the present paper is to extend such analytic descriptions 
of leading-order amplitude and phase fluctuations in the dispersion-managed 
coupled breathers derived from the Gaussian pulses as a solution of the cou- 
pled NLS with a perturbation term of a noise sourced by amplifler. We use an 
extended average variational method to reduce the complicated coupled NLS 
with a periodic dispersion map to several set of complicated nonlinear ordinary 
differential equations (ODEs), which accurately predict the completed coupled 
breathers dynamics. In order to flnd the Hamiltonian dynamics of this sys- 
tem, the ODEs arc derived by reducing and finding the parameters amplitudes 
rjj {Z) , and chirps Pj {Z) in which frequencies fj (Z) and center position of the 
pulses Tq. [Z) have been set to be identical with zero values. 

The paper is organized as follows. In Section 2 we introduce the appropriate 
extended coupled NLS with a perturbation term and its relevant physical 
parameters. In Section 3 we derive the Lagrangian of the nonintegrable coupled 
NLS with the perturbation terms using variational method, and find that in 
order to make the method work in this system, the parameters frequencies 
fj (Z) and center position of the pulses Tq . (Z) must be chosen to be very 
small (w 0). Section 4 transforms the extended coupled NLS equations from 
Lagrangian to Hamiltonian formulation for which the solution of the equations 
of the coupled breathers dynamics can be expressed in terms of a modifled 
quadrature. Section 5 contains of phase-plane dynamics of the derived coupled 
breathers dynamical system. Section 6 is the reduced results of section 5 
related to the single NLS provided by Kutz, et al., and the unperturbed coupled 
NLS. These system solutions arc found by reducing some parameters of the 
coupled NLS which has the perturbation term. A conclusion of the results is 
provided in Section 7. 

II. Formulation of the Nonintegrable Coupled NLS with 
a Perturbation Term 

A more complicated propagation of optical pulse in DM fiber is described 
by a non-integrable coupled nonlinear Schrodinger equation (NLS) with peri- 
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odically varying dispersion a {Z) and some additional perturbation terms: 



+ "^f ^ + HQil' + \Q2n Qi = iG {Z) Qi + ieiRi 
+ + (iQil' + \Q^V) Q2 = iG (Z) Q2 + ie2R2 

where Qi and Q2 are the electric- field envelopes normalized by the peak field 

22 M " 

powers |i?oi| and I-E02I , G (Z) = —T + Ga S {Z — nZa) represents the 

n=l 

eflFects of the fiber loss and the amplification, where F is the normalized fiber 
loss, Ga (Z) (= rZa) is the gain of amplifier, the following term^ 

ejRj = exp {i [fj {Z) (T - To, (Z)) + } , j = 1, 2 

are noise factors of the amplifiers (1 e) where raj, fj {Z) ,Toj {Z) and {Z) 
are complex amplitude of the noises, frequencies of the noises, center position of 
the pulses, and the phases of the pulses, respectively and we have set fj (Z) ,and 
(Z) to be very small, and aj (Z) is a pc;riodic dispersion map. To follow up, 
the previous work of Kutz, Holmes, et al., in their single NLS case, we define 
the variable T as the physical time normalized by To/1.76, where To is the full 
width at half maximum (FWHM) of the pulse (Tq = 50 ps). The variable 
Z represents the physical distance dcvided by the dispersion length Zq, which 
corresponds to the average value (£)) of the dispersion map a (Z). Thus we find 



XqD \1.7Q 



which gives 



for the peak field intensities of the one-soliton solution corresponding to the 
average dispersion. Here n2 = 2.6 x 10~^^ cm^/W is the nonlinear coefficient 
of the fiber, A^ff = 55 /xm^ is the effective cross-sectional area of the fiber, 
and Ao = 1.55 /im and c are the carrier's free-space wavelength and speed of 
light, respectively. In the present analysis, we involve the periodic intensity 
fluctuations that are due to fiber loss and amplification because they have very 
small perturbations to the average pulse dynamics. Thus the peak powers given 
in eq.(3) are almost the same as the path- average peak power of the propagating 
pulses. Note that we also consider long-length scale fluctuations in intensity that 
are due to gain-shape changes and saturation effects of the small noises in the 
amplifiers. 

The parameters that are of greatest interest correspond to the particular 
dispersion map that is being used. It is helpful therefore to consider the specific 
case of a piecewise constant step-function dispersion map such that 
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ai (Z) = <72 (Z) = ^{ D+, 



< Z < 
1^ < z < i— + ^ 

2 Zq 2 Zq Zq 



(4) 



where a (Z) ^ a {Z + P) , D+ > > D_ are the dispersions in each segment of 
fiber (in picoseconds per kilometer. nanometer), Z± is the length of each fiber 
segments (in kilometers), and the average dispersion is simply given by 



D = 



D+Z++ D_Z. 
Z-+Z+ 



(5) 



We also follow the values of the parameters {D = 0.2 ps/km.nm, Z_ = 450 km, 
= 60 km, D- = —2.1 ps/km.nm, and = 17.45 ps/km.nm) which are 

used by Kutz et al., in their explanation related to the single NLS. 
In order to simplify eq.(l), we use a transformation as follows 



Qj (Z, T) = VnZ)Uj (Z, T), j = 1, 2. 



We then find 



dU2 



where the reduced noise factor used in this investigation is 



T{Z) 



Uif + \u2n U2 



y/r{z) 

if-2R2 



(6) 



(7) 



y/Vj^aj exp 



in which the amplitude factors ,yf]~ used here are to reduce the influence of 
the noise in the chirp dynamics and T{Z) is chosen by following the work of 
Okamawari, et al.^ in the case of their single NLS of eq.(7). 



r{Z) = T(0) exp 



! ^ G{Z')dZ' 



(8) 



and the average values of T(Z) and a {Z) are assumed to be unity in the nor- 
malized unit. The imaginary term in the perturbation term can also be made 
as a real number by choosing, for instance, the complex amplitude of the noises 

raj = exp{-0.5i(l)j (Z)). 



III. Variational Method of the Coupled NLS 

The fundamental physical principle involved in this approach is derived from 
Lagrangian methods or Euler-Lagrangian equation for classical mechanics writ- 
ten by Goldstein, H.^ This requires the dynamics to minimize a certain energy 
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integral closely related to in Hamilton's principle. We can restate the extended 
coupled NLS eq.(7) in terms of its variational form by defining the Lagrangian 

oo 

£{UjU;)dT, J ^1,2 (9) 

— CSO 

where £ (UjUj) is the Lagrangian density for the eq.(7) given by 



£{U,U*] 



2 

E 



, du* 



dz 



^ dz 



a{Z) 



dT 



T{Z)\U^ 



(10) 



In order to reproduce eq.(7) and its complex conjugate, the Lagrangian density 
must be determined by requiring the following equation 



5£{UjU*) S£{UjU* 



SU* 



5Uj 



= 0. 



(11) 



In the reduced model developed here, we assume that the fundamental dy- 
namics occurs owing to changes in amplitudes and quadratic phase chirps and 
they depend on very small noises sourced by amplifiers. Further, we assume a 
solution of eq.(7) that has an extended Gaussian shape as follows 

Uj {Z, T) = A,&,VT7-exp {- {k^t^^tY + i [fi^K^T^ + ^£1 j } , j = i, 2 

(12) 

where rj^ {Z) and Pj {Z) are free parameters that correspond to the amplitude 
and quadratic-chirp dynamics, respectively. We then find the total energy pulses 
as follows 



oo 

J (\Q,{Z,T)\'' + \Q2{Z,T)\^)dT 

— OO 

= j [{hi exp (-2K?r2)) + [bl exp {-2kIt'^))] dT 



(13) 

Kl V 2 K2 V 2 

As with the governing eq.(l), this ansatz conserves the total energy of the 
pulse, i.e., 7 (iQil' + IQ2I') dT = (^^) + (iM^) = constant. More 
generally, any localized ansatz for which the amplitude and width are related 
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as in eq.(12) can be utilized. Related to the case in the same width (FWHM) 
as the one soliton solution, we find + ^ 

In general, approximating the eigenmode of eq.(l) with the dispersion map 
a (Z) by a Gaussian profile leads to errors in the reduced dynamics. However 
these errors can be reduced by an appropriate choice of the parameter Aj. To 
extend the previous work of the single NLS provided by Smith et al.^, we define 
the modified amplitudes-enhancement factor Aj as follows 




J = 1,2 (14) 



The amplitudes-enhancement factor aj ^ 0, if we take an arbitrary ^ or 
(Aj = ^ 1 in order to characterize a more complicated the Gaussian pulses 
dynamics. This problem can also be easily discussed using a numerical compu- 
tation and simulation. However, in order to simplify our analytical calculation 
we choose Ai= A2 = A and then find ki = k2 = ttI if h\ = h\'^ 1. 070 2. 

After substituting the pulses ansatz of eq.(12) in the Lagrangian equations 
of eqs. (9) and (10) and evaluating the Lagrangian, we then find 



dcj)^ jZ) _ r{Z)^r]j 
dZ 2 



^CTj (Z) 



-2T {Z)r^,V2AlAlblhl^l^ (exp [{nl^l + ^r^l)]) 



] ^20F 
where we have chosen 



01 ( ' 



or 



and found 



Re ^(kj-tj^)^^ = and Im ^(KjTy^)^^ > 0, 



(15) 



"^■ = E^K-^2,), (16) 



where 
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-u,2,)= exp [-^j [exp {i<f>j {Z)) - l] (17) 

Applying the variational method with the Gaussian trial function of eq.(12), 
we obtain that the related parameters of Xo^ ,aiid fi in the noise terms are very 
small and equal to zero. 

Although the complicated variations with respect to the free parameters T]j , 
and /3j can also yield the appropriate ordinary differential equations (ODEs), 
we will extend the Lagrangian formulation of the perturbative coupled NLS 
and will then try to develop the Hamiltonian formulation of these equations by 
means of a Legendre transformation. This method reveals a conserved quantity, 
which is utilized in the solution of the nonintcgrablc coupled NLS equation. 
The validity of the variational method and its associated dynamics are also in 
order related to the certain Gaussian pulse ansatz such as provided in eq.(12). 
Then, we can produce the correct pulses evolution. This certainly is not true in 
the case for which self-phases modulation dominates dispersion and temporal 
side lobes are generated, i.e., if Aj are too large. 

Although Kutz et al. had already said that the variational method was 
incapable of describing or capturing the slight growth of continuous spectra, 
which is generated by the dispersion-management perturbations. We improve 
that it is possible if wc can reduce all of the additional parameters such as the 
frequencies fj (Z) and center position of the pulses Tq^. (Z) to be identical with 
zero. 



IV. The Hamiltonian Dynamics of The Coupled NLS 

There are various ways and manners in deriving the solution and explanation 

of the amplitudes and chirps of the nonintegrable coupled NLS. However for 
those cases in the nonintegrable systems, almost all of the methods cannot 
work to derive the solutions. So, the only useful method we are going to use 
is the Hamiltonian formulation of the dynamics. The Hamiltonian system for 
the system is useful since solutions can be derived and represented in terms 
of a modified quadrature. The Lagrangian system of the preceding section is 
brought to Hamiltonian from via a Lengendre transformation. We can directly 
define the related canonical coordinates (3^ with its conjugate variables as 
follows 

dL ^ 

P^^ = T77irY = Z^^ i = l'2 (18) 



d 



which represents the canonical momentums. The Legendre transformation de- 
fines the Hamiltonian 

h{p„p,^)=p,^^-l(p„p,^), (19) 
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where we must change the Lagrangian in eq.(15) into the Lagrangian in the new 
canonical variables and p/}. and then substitute it into eq.(19), we obtain 

.7 = 1 V ' / 



+2^1^ (Z) (p^ J (p^J vlvlblbl exp {kWi + i4vl) 
( i20r 



[ra.^i (/3i,P/3j +ra,r!2 (/?2,P/3j} , (20) 



where Q.j {jij,Pi3^ are defined as follows 



The Hamiltonian can then be calculated with eq.(19), and the evolution equa- 
tions for r]j dan (3^ are determined from the following Hamilton's equations 

dZ ~ dp0. ' dZ ~ d/3j ■ ^ ' 

In terms of r]j , we can evaluate the evolution c;quations in eq.(22) to be four 
ODEs stated the dynamical system of the Gaussian pulses as follows 



^ = 2a,iZ)^l[rjt-Pl]-^riZ)4Albl 

-2T {Z) rj^rjlAlblm (l + 2rilKl) exp (k???? + kItjI) 

Ki \ i2\/2rai (9 

— i ll 



(23a) 



dp2 

dz 



= 2a2{Z)4[vi-Pl]-^vlAiy-, 



-2T {Z)r]^r]lAlblK2 



27j14) 



exp (nl-ql 



2 2\ 



(23b) 



and 
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where 



= (eK(^)) - 1) e ^ (24a) 



and 



90^ ' 27^2 



(e('^.(^)) - 1) e "^^ (24b) 



The results of the ODEs in eqs.(23) describe that the coupled pulses dy- 
namics depends on the chirps interactions contributed to the pulses amplitudes. 
Loss and gain fluctuations related to the center position of the; pulses and the 
frequencies of noises are set very small to this formalism. Hence, the system 
of equations decribes the nonlinear enhancement amplitudes and width fluctua- 
tions along with the modifled quadratic-chirp variations as the pulses propagate 

through a given dispersion map. Note that the resulting ODEs system is iden- 
tical to that derived by taking the variations 

^=0 and § = 0. (25) 

Since the resulting nonlinear ODEs system of eqs.(23) is the Hamiltonian, 
and aj {Z) is piecewise constant, we can utilize the fact that 

and construct level sets of a scaled version of H rewritten in terms of rjj and 
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C = 



Kl 



1 

nl dz 



^9 ^ 



( .2''i 

1 dZ 



A2_L 



> - fTl (Z) 



1 dp^ 



/3? 



iflr^2 (^) K2 



. ( r{z) 



■) M 



+ 



+2T (Z) A2Ki7yir?26i (exp7(Kf r?? + Klrjfj] ) 



In order to reduce the eq.(27), we choose 



2k2 



and 6i 



■(4-2^i) 



Kl 



2 / \ ^iVT 

The eq.(27) can be more simpUfied by choosing 



b^. = = b^ K 1.070 2, 

and [ki = K,2 = K = -j^) 

Al = Al = l. 

We then find the general solution of (3i (r/i, ??2) P2) follows 



(27) 



(28a) 



(28b) 



(28c) 



= ±771 {71 [Vi + V2] - 72 - 73 {vl + f ) + 74^i%e['''(''^+''^)] + 75 " } % 

(29) 

where 



7i 

72 
73 
74 



2T{Z) 
a {Z) K\pK 

( Q I rf0l(g) d^2(£l\ 

1 (for ai (Z) = (72 (Z) = a) . 

4V2T (Z) 

C7(Z)V^' 



(30a) 

(30b) 

(30c) 
(30d) 
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and 



75 = (^ ^^z)l\/T{Z) ) ^""^"^ +^a.^2 (30e) 

where fli (/3i,P/3^) and 02 defined in eq.(21). 

V. Phases-Plane Dynamics of The Perturbed Coupled NLS 

It is often convenient and insightful to plot the phase-plane dynamics of a 
planar dynamical system such as eqs.(23). This geometric representation of the 
solution treats the distance Z as a parameter while plotting the amplitudes 77^ 
and chirps (3j variables on the x and y axes, respectively. It elucidates some 
fundamental aspects of the coupled breathers dynamics, which will be clearly 
shown in what follows. 

The phases-plane dynamics are markedly different depending on the sign 
of the dispersion. Thus we consider separately the anomalous- and normal- 
dispersion regimes. Before seeking the fixed points of eqs.(23), we will simplify 
the very complicated equations in order to get the phases-plane dynamics con- 
sisted of (3i {r]i), Pi (772), 02 iVi) and (^2 iVi) ■ This will easily work by choosing 

I32=pAZ)Pi, (31a) 

and 

V2=P2{Z)Vi, (31b) 

or by setting (Z) = p2 (Z), we suggest and make our solution in the simplest 
forms: 

p(Z) = ^ = ^= const. (31c) 

Pi Vi 

After substituting eq.(31c) into eqs.(23) and eq.(29), we then find 
Pi = ±Vi l^laVi - ^ - IbVl + Icnl exp [kWi (1 + p')] + ^ - ^} \ (32a) 

P. = ±.2|^.2-^-^.^ + %^(e"^^^"^>0+^-^4'' 

[ P"* P^T p^ p^ \ J p'^T p^T ) 

(32b) 
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/32 = 



2 7c^i(e'='('+^)''^ 

, „ ; "faV2 72 76^2 , V / , 75 V2 

-^'12 \ 2 ' 2 ' 2~ 



(32c) 



/32 = ±77i <; p^7„??i 



T r T 



(32d) 



where we have reduced the eq.(32) by choosing r as follows (cti (^) = cr2 (^)) 



r = 2, 

and the parameters 7^, 7(,,7g,and the revised 75 are 



(33) 



7a 
76 

7c 



7i [1+P(^)] 
2 

2 ' 

74P(^) 



(34a) 
(34b) 
(34c) 



and 



75 = 



\a{Z)t^^^/Y{Z) 



I [fii + O2] 

^(#1(2)+) ^g(i</.2(Z)) _2' 



, for = = r^, 

(34d) 

or by choosing = exp [— 0.5i(/)j (Z)] in the eq.(34d), we can reduce the 
imaginary and then find 



75 



[sin(0.5(?ii) +sin {O.bcl)^)] 



(34e) 



where ( S {Z) > 1) . In eq.(34d), we have suggested the following relationships 



y/nvab'^ (e'«^i + e'«^2 - 2) 



1 6^ e'*"! + 6**^2 _ 2 \ 
LambertW ( -i^/^^ra-y= — (L)j , 

(35a) 
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where 



and 



and the rule appeared in eqs.(35a) and (35b) must be used to our entire solutions 
in this section. If we choose S (Z) = 1 or the chirp /3j « , this then means that 
the pulse width is minimum at the center of unit cell. However, the mean-square 
frequency shift is not minimum. We can also suppose that a noise (in 75 term) 
sourced by the amplifier which has a small frequency causes a small change in 
the coupled breathers dynamics of the Gaussian pulses. However, we still can 
reduce the noise by enhancing the amplitude of the pulses. 

The aim of the use of a relationship in eq.(31c) is to make the phases- plane 
dynamics of the coupled breathers are in stable conditions. A first order ODE 
in terms of TJialone can then be derived by substituting eq. (32a) into eq.(23c): 



^[ 1^T{Z) ) 

X jla^/i - ^ - + exp (1 + P^)] + ^ - ^} '(36) 

Upon integrating this equation with initial conditions r]y = r]Q^ at Z = 0, we 
find 



7 T {laV, - ^ - 7,vi + l.vl exp [nW, (1 + pI)] + - ^}"^ dr^. 



(37) 



Another first order ODE in terms of alone will easily be derived using the 
same manner. Although we already get an analytic solution for the pulse dynam- 
ics in eq.(37), the integration is still very complicated due to some closed-form 
expressions. The simplest way left to explain the coupled breathers dynamics 
is a geometrical interpretation of the ODEs in eqs.(23). 
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Finally, we come to obtain the fixed points of eqs.(23) by setting the right- 
hand sides equal to zero. There are three fixed points which can be proved 
by figuring the phases-plane dynamics consisted of p-^ {r]-i), Pi (r?2), P2 {'Hi) and 

I. (3. = 0, r?,- = 0, (38a) 



II. f3j = 0, and 

an an - 0,41 exp (t?i) 031 - 071 exp (t?i) 031 



III. Pj = 0, and 

_ C21C71 

C7icn - C41 exp C3ii?i - C71 exp (i?!) C31 

where i?! is a root of 



—2e® 0710110410316^ — O21O71 + ©a^iO^i}, (39a) 

and i?'i is a root of 

{e^^Oc^iC^i + 2e2®c7iciiC4ie2 + e^^c\yclyQ'' - 2e'^e^^iCnC3i 
-2e®C7iCiiC4iC3ie2 - ciic^i + ec^icfi}, (39b) 

and the parameters r?oi, '?02i an, 021, 031, 041, 071, cn , C21 , C31 , C41 ,and C71 (for 
Ai = A2 = A and ai (Z) = a2 (Z) = a {Z)) are 

%1 ~ 2k1(1+p2)o--2k2 exp(i?i)[2T(Z)pA262K]^)i-K2(i_|_p2) exp(i3i)[2T(Z)pA2&2„] , 

„ _ 

2k*(i+^)<t-2k2 exp(iJ'J[2T(Z)pA2b2K]tf'^_«;2(l+^) exp(iJ'J[2T(Z)pA2b2^] ' 



2f7 (Z) K^, 

A^T{Z) fe2 



Oil 


= Cn 


a2i 


= C21 


asi 


= C31 


041 


= C41 


071 


= k2 


C71 


= 



V2 ' 
^'7^ „/i2i,2,. 



„2\ 
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A linear-stability calculation about each of the fixed points gives insight into 
the general coupled breathers dynamics of the system. In particular, the fixed 
points (77^, /3j) = (0,0) are degenerate so that their linear stabilities cannot 
be determined (sec, for examples, figures 1, 3, 5, and 7). However, the linear 
stability of the fixed point II and III can be determined by studying the pulses- 
plane dymanics of the ODEs in eq.(23). Linearizing about these points by 
setting 



(40a) 



f3j = + f3j 



(40b) 



yields the linear system and it can be rewritten in a matrix form as follows 



where 



/ dvi \ 

dZ. 

^ dz ' 



f 








Mi3 









\ 













M24 








M31 















Pi 







M42 








) 




\ P2 J 



(41) 



Mi3 = 



-CTl {Z) K?7oi 



.(#i(Z)) 



1) Kb^E{Z)rjoi 



(42a) 



24 



-6 ^/-£72 {Z)K'nQ2 



ira^ V'' 



1) KlPZ (Z) ,i,y. 



2^/T{Z) 



(42b) 



M31 = 8a,{Z)K'vli~^r{Z)AlbWoi 

-8T (Z) pAlb^Krjl, (1 + 27jI,^) e('='(i+^')''oi) 

-8T (Z) pA262^3^5^g(K^(l+p^)0 

-4T (Z) pAlb^K^^rf^^ (1 + 2if^^K') (l + p^) e^^'i^+o'W (42c) 
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M42 = 8a2(Z)«\^2-^T^i?'''?02 




The eigenvalues of this system (eq.(41)) determine the stabihty near each of the 
fixed points and are given by 



The eqs.(43) show that the fixed points II and III ((rji.Pi = 0), (?72,/3i = 0), 
{r/i, ~ 0) and (772, ^2 = 0)) are the centers if the sign of the dispersion (a) is 
positive. We will then investigate both the anomalous- and normal-dispersion 
fibers in what follows. 

A. Normal Dispersion: (ai = C72 = cr) = < 

We start by considering the phases-plane dynamics in the normal dispersion 
regime. In this case, we can rewrite a = — ^-^^^ < 0. We then consider the 
location and stability of the three fixed points of eqs.(23) shown in eqs.(38). 
Fixed point I has already been determined to be generate, and fixed points II 
and III lay in the left half-plane since iq^. < 0. Because we consider only values 
of r]j > 0, only the fixed points at the origin are relevant, the shapes of the 
plane are the same as homoclinic orbits, which emanate and terminate in the 
origin. 

B. Anomalous Dispersion: {ai — a2 — (j) = ^ > 

In the anomalous-dispersion regime, the phases-plane dynamics are signifi- 
cantly different than that of the normal-dispersion regime since fixed points II 
and III , given by eqs.(38b) and (38c) are located at 




±V(M3iMi3), 



(43a) 
(43b) 



±V/(M42M24) 



Pj =0, and 



II. 



>o, 



D 



\kP- exp(i?i)[2T(Z)^jA252re]i?i-K2(i+p2)gxp(i?i)[2T(Z)pA262K] 



(44a) 



and 



16 



f3j =0, and 



III. 



1(1+ i)d, 

^ 2k2 cxp(-d[)[2T{Z)pA^b^K]-&[-K^{^l + -^) exp(i9;)[2T(Z)pA2b2K] 



>o, 



with the four eigenvalues ^for a = 



(44b) 



Al,2 
A3, 4 



±V(M3iMi3), 



(45a) 
(45b) 



In this case, the phases flow for rjj > are partially determined by the 
three critical points I, II and III . It can be more obviously understood if we 

replace C = - (^^^^ + ^^^^) in the eqs.(32) and = exp (-O.Sic/-^. {Z)) 

in eq.(30e), We then give rise to the separatix 



01 = ±Vi \ laVi + IcVl exp [k't?? (1 + p')] +Y-[^b+^]vl\ , (46a) 



Pi = 



1 + 276 ^ 2 
2p4 ) ^2 



(46b) 



02 = ±'72 < 



2 / K^(l+^)vl 
laV2 , "'V y , 75 ^1 + 27m 2 



2p2 



'?2 



(46c) 



and 



= ±r]i{p'laVi+p'lcVl{e'^^'^^^'+^^^) + 



(46d) 

which have a cusp at the origin (/3^. = r/j. = 0) . The solut ions outside these 
separatix eventually flow into the origin (homoclinic orbits), while those inside 



17 



the scparatix arc periodic. In cqs.(46), there is an indicator 75 represented the 
contribution of a perturbation term caused by a noise sourced by ampUfiers 
in optical fiber and it can influence the amplitudes and chirps of the pulses 
dynamics. The geometrical figures and its analysis of this system of cqs.(23) 
are not provided in this paper due to a complicated interpretation. We will 
then show the simplest explanation of the unperturbed dynamical system in 
the following section. 

VI. The Hamiltonian Dynamics of The Unperturbed Cou- 
pled NLS 



The solution of (3-^ in eq.(29) can be reduced to that in the single NLS 

which has no the perturbation terms provided by Kutz, et al. by setting = 
f]i = 4>i = ^1 (■^) = 1 and tjRj = 0. On the other hand, we conclude that 
the terms of the Lagrangian, the ODEs and the (3 solutions of the unperturbed 
coupled NLS : 



; dUi 
' dZ 

:dU2 



2 dT^ 

<y{z) d^u2 
2 'SW 



\Ui\ 
\Ui\ 



"az ' ^5^7Fr2- I I \U2\ ]U2 

reduced from the results in section 5 are, respectively, as follows 







(47) 



= E 



I a/3, 



2772 dZ 



+ 




Gj (Z) K^j 



dZ 




-2^ViV2AlAlblbl (exp [{kWi + 44)] ) , (48a) 



^ = 2a,iZ)nl[4-Pl]-^rjlAlbl 

-2r]2vlAlblKi (1 + 244) exp + 44) 



(48b) 



dZ 



2a2{Z)4 [vt-4]-^4Albl 
-2rii4AlblK2 (1 + 244) exp {44 + 44) 



dr] 



1,2 



(48c) 
(48d) 
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and 



(49) 

The simplest stable solution of the equations for the amplitudes r]j and chirps 
Pj derived from eq.(49) will be found by choosing 7/2 = PVij and $2 = P/^i) 



/3i = ± 



V2 

pV2 
{72) 



?eK(i+''^K]-ry2}\ (50a) 



2 

7i [1 + P] — - 72 - »?2 + 74— e 



-"^2 



9 

'72 



, (50b) 



/J2 = ± ( ^ ) {71 [1 + p] r?i - 72 - + l.pnle^'^'^'^o^^'^ - r??} % 

(50c) 



and 



/32 = ±^^7i[l + p]^-72-r?i + 74fe 



— -^2 



(50d) 



where (for cri (Z) = (T2 (^) = cf [Z) and 61 = 62 = h) 

2 



7i 

72 
74 



a (Z) Ky/n 



' dZ dZ 



a (Z) k2 



4V2 



a(Z)0F' 



The substitution of eqs.(50a) and (50d) into eq.(48d) yields the 
first ODEs in terms of r/^ and 772 , respectively: 



and 



(51a) 

(51b) 

(51c) 
following 



(52a) 



dV2 _ -j- b ^(t{Z)k7J2 
dZ ^ 2 



7i [1 + P] ^ - 72 - + 2f^e 



Ik 

p2 



(52b) 
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Upon integrating these equations with initial conditions tlj = rjQ. at Z = 0, we 
find 



T2 



/ 



h^y/^a {Z) ktjI 



(53a) 



or 



T2 



V2 

I 

V02 



7i 







^ 5 — ^ 

P 


^1 



f^y/^a {Z) nrii 



(53b) 



where a (Z) is a fixed value, i.e., cti = (T2 = or depending on whether the 
coupled pulses arc propagating in the normal or anomalous regime, respectively. 
The integrations in eq.(53) are very complicated expressions, which must be 
transformed to give t]^ = rj^ (Z). However, we will explain it elsewhere. The 
only simplest manner in describing the coupled breathers dynamics is by using 
a geometrical interpretation of eqs.( 48b), (48c) and (48d). 

The fixed points derived from a perturbative coupled NLS are usually very 
complicated in terms of the roots. However, in order to explain the phases-plane 
dynamics of the coupled breathers, we come to obtain the unperturbed fixed 
points of eqs.(48) by setting the right-hand sides equal to zero. Here, we find 
three fixed points reduced from that in eqs.(44): 



I. /3, = 0, 



T], = 0, 



(54a) 



II. /3, 



071011 - 041 exp {'di) asiXfi - an exp (t?i) asi 



III. Pj = 



_ C21C71 . . ^ 

C71C11 - C41 exp (i?i) C3ii?i - C71 exp (i?J C31 



where i?! is a root of 
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— 2e®a7ioiia4ia3i0^ — a2i'^7i + ^O'nO'ii}^ (55a) 



and 'd'l is a root of 



„2e 



26 „2 2 



-26*^0710110410316^ - O21O71 + 6O71O11}, 



(55b) 



in which the roots i^iand i}[ arc, respectively, {)i = 1. 316 5, and t)'i = 6. 189+3. 
136i i{A = l,a = 1, and p = 0.03, and the parameters r?oi, ryo2, «iij 0215^31, 041, 
071, cii, 021,031,041, and C71 (^i = A2 = A) are 



_ ^A2fc=K=(l+p2) .gg 

%1 ~ 2k*(1+p^)ct-2k^ exp(i3i)[2pA^6'-^K]i?i-K;^(l+^)^)exp(i?i)[2pA^6^K] ' ^ 



2k4 ( 1+ ^ ) o-_ 2k2 exp(i?'J [2pA2 62 «] - k2 ^ 1+ _^ ) exp (iJ'J [2pA2 62 ^] 



(56b) 



an 


= Cii = 


2a {Z) 


a2i 


= C21 = 


A^^ 
^/2 ' 


^31 


= C31 = 


2pA'^b'^K, 


^41 


= C41 = 




071 


= ^{1 


+ p'), 


C71 







Linearizing of eqs.(48b), (48c) and (48d) by setting t]j = rj^^ + rjj and Pj = 

0+(3j, we find the hnear system by replacing ejRj = or rm = 0, and T (Z) = 1 
in eqs.(42) and it can be rewritten in a matrix form as follows 



where 



/ ^ \ 

dZ 
Ml 

V Ma / 











Mi, 









/ Vi 


\ 













MU 












A^31 















Pi 




\ 





M12 








) 




\ /?2 


/ 



= -62 



cr {Z) K'no^ 



(57) 



(58a) 
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MU = -b^^c7{Z)KVo2 (58b) 

-8pAH\4, (1 + 2772,^2) ei'^'i^+P^yoi) 



(58ci) 



The eigenvalues of the system (eq.(57)) determine the stabihty near each of the 
fixed points and are given by 



Al,2=±{ 



and 



V2'^ '01 ' 8»JoiKe''iptfi+4»7oiKe''ip+4r)oiKe''ip3 + y2+V2p2 



+ 



32 



per 



(1+P=) 



(8r7oiKe''ip!?i+4r;oiKe''ip+4»7oiKe''ip3 + \/2+\/2p2) 



'/ill 



(1 + /'^) 



p3 + V^+V^p2) 



32 



+■ 



(8r;o 1 Ke'' 1 pi?i +4t7o 1 fte'' 1 p+4j7(,i Ke'' 1 p3 + V2+ V^p^ ) 
16^pa=6^«^Ji(l+2„g,K=)(l+p=)e('''(^+''')''oi)l (l+p^) ^ 



(87701 «e''i pi?i +47701 ree-'i p+4r)oi «e''i p3+\/2+\/2p2 ) 



(59a) 
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r 12 „2 



32 



262KSo2(l+2';o2'*^)e 



,2^ (^=(^+;^).g2) 



(l+p=) 



p(8r)02 Ke''i p!?i +4r)o2 «e''i p+4r)o2 Ke**! p3 + \/2+\/2p2 ) 



+ 



+ ■ 



p(8^ 



, Ke''i pj?i +4r)o2 Ke''i p+4?)o2«e''i p3+\/2+\/2p2 j 



16^a=6^«^j2(l+2'7^2«') (1+^) e("' "'^ ) ""^^ 



p(8??02KeOi pi>i+4??„2K:e*i p+4r;o2Ke*'i p^+s/2+^/2p^) 

(59b) 

The eqs.(59) show that the fixed points II and III are the centers, regard- 
less of the sign of the dispersion. We will then investigate both the anomalous- 
and normal-dispersion fibers in what follows. 

A. Normal Dispersion: a = S=- < 

In this case, we can rewrite cr = — ^-^^ < 0. We then consider the location 
and stability of the three fixed points of cqs.(48b), (48c), and (48d) shown in 
eqs.(54). Fixed point I has already been determined to be degenerate, and fixed 
points II and III lay in the left half-plane since rj^. < 0. Because we consider 
only values of 77^ > 0. only the fixed points at the origin arc relevant. The 
phases- plane dynamics are shown in figures 1, 3, 5, and 7 along with the arrows 
in phases superposition space. We find that the shapes of the plane are the 
same as homoclinic orbits, which emanate and terminate in the origin. 

B. Anomalous Dispersion: a = = > 

In the anomalous-dispersion regime, the phases-plane dynamics is signifi- 
cantly different than that of the normal-dispersion regime since fixed points II 
and III , given by eqs.(54b) and (54c) are located at 



II. 



/3,. = 0, 



2K4(l+p2) jL_2K2exp(l?i)[2pA262K],?l-K2(l_|_p2)exp(??i)[2pA262K] 



>o. 



(60a) 



and 



/3, = 0, 



III. 



2k;4(i+-^^ iii-2K;2 eyi-p(^i)[2pA'^b'^ k]'&i-k'^ exp(i?i)[2pA262«- 



>o, 



(60b) 
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with the four of eigenvalues 



Ai,2 = ±{ -8(^)%^v^- 



.3 4 
Vol 



— (E+Y 



+■ 



8r;o 1 Ke'' 1 pi5 1 +4 jjo , Ke'' 1 p+4?7(, 1 KC.'' 1 cr' + \/2+ V2p^ 



+ 



'32V¥p(^)%^«%J.e(^^(^+''^)''o01(l+p^) 



01 Ke**! pi?i +47701 Ke"! p+4r;oi 



ip3+\/2+\/2p2) 



+ ■ 



(8))oiKe*'ipi9i+4»7oiKe''ip+4r;oiKe''ip3+v^+72p2) 
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p(^)^^K^,^i(l+2r,i5iK^)(l+p^)e('''(l+-')''oi) 



(8t7oiKe''ipiSi+4i7oiKe''ip+4?7oiKe''ip3 + V2+V2p2) 



(61a) 



and 

•^3,4 



V2 



+ ■ 



8jJo2Ke''i p'i5i+4jJo2Ke''l p+4r;(,2Ke''l p3+V2+\/2p2 



p(8)7o2 Ke'' 1 pi?i +4j7o2 kb^ i p+4j7o2 Ke'' i p3 + ^2+ V2p^ ) 



+ 



32yf(^)\=«So%e(^^(-^)''o 



p(8_))o2 1 piJi +4)7o2 fee'' 1 p+4))02 fte'' 1 p3 + y2+y2p2 ) 



32Vf(i^)\^«SLe(~^(-;^)''-) 



+ ■ 



p(8))02fte''ipi?i+4))02Ke''ip+4))02Ke''ip'Hy2+y2p2) 
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yf(^)\^K%J,(l+2r4K^)(l+-^)ei 



p3 + V2+V2p2) 



(61b) 

The phases-plane dymamics arc shown in figures 2, 4, 6, and 8. In this case, 
the phases flow for rjj > are partially determined by the three critical points 
I, II and III . It can be more obviously understood if we replace 



C = - 



dZ dZ 
in the eqs.(50), We then give rise to the separatix 



- -i-Jh 



/3i = ± 



;+2]gr,2 

pV2 "I D+Kv^ 



D4V2', 
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+ 



exp[K^(^ + l)7;^] 



(62a) 
(62b) 
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(62c) 

and 

which have a cusp at the origin = rij = O) . The solutions outside these 
separatix eventually flow into the origin (homoclinic orbits), while those inside 
the separatix are periodic. 

From figures 2, 4, 6, and 8 in which the parameter p has been set to be, for 
instance, p — 0.03, wc then find exactly that the fixed points of (rii,Pi — 0), 
(^2)/^i = 0)) (^i)/52 = 0) S'lid {r]2,P2 = 0) ^■t the center of their graphics as 
follows 

(?7i,(/3i=0)) = ( r,i,(/32= 0)) = ( 0.75 , ), (63a) 

and 

( V2, iPi = 0)) = ( r?2, {P2 = 0)) = ( 0.0225 , ) . (63b) 

The points and the phases-plane dynamics of the ODEs eqs.(48b), (48c) and 

(48d) improve that our suggestion in eq.(31c), p{Z) = ^ = ^ = const., is 
actually proved. This means the coupled breathers dynamics depent on the 
certain relationship of the chirps and amplitudes. The choice of an arbitrary 
p{Z) makes the certain fixed points of ('7j,/3j) obeyed by eqs.(54). These fixed 

points are found in + '^''^^^^^ + ^^g^^ ~ 0.605. In the physical meaning, the 

interactions between both chirps depend on a linear constant p{Z). However, 
the magnitude of the Gaussian pulses amplitudes is not only directly influenced 
by the chirp itself but also by the interaction of both chrips. 

On the other hand, we, for instance, can set the values of p (Z) and then 
find the the amplitudes by using a numerical method from their phases-plane 
dynamics graphics. 

A geometrical representation of the solution and its associated dynamics 
behavior of eqs.(48b), (48c) and (48d) appeared in figures 1-8 show the super- 
position of the flow of phases ((^^ and (/ij). According to the figures, the con- 
struction of a periodic dispersion-managed coupled breathers solution from the 
dispersion map depends on the critical values of Aj. In this periodic solutions, 
the pulses first evolve according to the first segments of the dispersion map for 
which cTj < (figures 1, 3, 5, and 7). The dynamics then reverse flow according 
to the dispersion switches sign as shown in figures 2, 4, 6, and 8, and then fi- 
nally are again subjects to the dynamics. This case is in a reflection symmetry 
of /3 • = derived from the choice of dispersion map and the resulting invariance 
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of cqs.(48b), (48c) and (48d) under the following reversibility transformation 
(rij,Pj,Z) (r]j,—Pj,—Z) . The qualitative desriptions are powerful enough 
when the exact values of the fixed points {VjTPj) (eqs.(63)) are directly found 
in the phases-plane dynamics graphics, and have been possibly proved in this 
paper. 

For the Single NLS investigated before by Kutz, ct al., all explanations about 
the Hamiltonian dynamics of the dispersion-managed breathers can be derived 
from the reduced ODEs of the coupled NLS in section 6. For instance, the 
two fixed points and the eigenvalues which are easily found by reducing eqs.(48) 
and (59) are as follows 

I. /3i=0, 7,1=0, (64a) 

II. /?! =0,V2 = 0, Vi = = Vio, (64b) 

anciii 

where f 6^ = -2^and p = 



and 

Ai^2 = H^aK^rjli (66) 



VII. Conclusion 



In conclusion, we have presented the recent results related to Hamiltonian 
Dynamics of Dispersion-Managed Coupled Breathers in Optical Transmission 
System. From mathematical viewpoints, the use of a variational method and 
its corresponding Hamiltonian formulation of the coupled NLS, the pulse dy- 
namics in a dispersion-managed optical system can be reduced to four nonlinear 
ODEs. By means of the conserved quantity of the Hamiltonian system, an exact 
solution of both (3j and rij can also be found in terms of a modified quadra- 
ture representations. The results of the coupled NLS with the perturbation 
term can also be reduced to both that system which has no the perturbation 
term and the single NLS case. From the physical viewpoint, we conclude that 
a noise sourced by the amplifier which has a small frequency causes a small 
change in the coupled breathers dynamics of the Gaussian pulses. We then get 
some phases-plane dynamics of the reduced system in several fixed points. We 
have also shown that the linearization and use of the Hamiltonian integral of 
the four ODEs describe the dynamics of the phases-plane and further suggests 
the possibility of constructing a periodic solution that depends on the initial 
amplitude-enhancement factor of the pulse power. Finally, we then find the 
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proof from the phases-plane dynamics that the coupled breathers dynamics de- 
pent on the certain relationship (eq.(31c)) and the parameter Aj, of the chirps 
and amplitudes. 

Thus, based on our results, there are some important open problems related 
to the another coupled NLS with a complicated loss and gain terms. However, 
the dispersion-managed vector breathers dynamics derived from the Gaussian 
pulses as a solution of the nonintegrable vector NLS : 



is in progress. As a result, the advantages for the application of the Gaussian 
pulses in optical communication system instead of the initial small chirps and 
pulses energy are still open. 
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